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Abstract 



Unitarity relates the total cross section for neutrino-nucleon scattering to the 
neutrino-nucleon forward scattering amplitude. Assuming the validity of the per- 
turbative expansion of the forward amplitude in the weak coupling constant, we 
derive a unitarity bound on the inelastic cross section. The inelastic cross section 
saturates this bound at a typical neutrino energy ~ 10® GeV. This imphes 
that calculations of the inelastic cross section that use current parton distribu- 
tion functions and lowest order weak perturbation theory are unreliable above 
this energy. 
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I. INTRODUCTION 



The phenomenology of ultrahigh-energy (UHE) neutrinos and their detection depends on 
the neutrino-nucleon total cross section, which has been calculated in the standard model in 
Refs. A striking feature of all of these predictions is the continued power-law-like growth 

of the cross sections with at the highest energies. This rise with Ei, is directly related to 
the very-low-x behaviour of the nucleon parton distribution functions (PDF's), which must be 
extrapolated below the regime of the current HERA data [|^ . It has even been argued p|-pT| that 



cosmic ray events will soon require a cross section even larger than that given by the standard 
calculations. 

On the other hand, unitarity relates the total scattering cross section to the forward scatter- 
ing amplitude. For neutrino-nucleon scattering, in contrast to electron-nucleon scattering, the 
forward amplitude can be determined and used to bound the total cross section. This bound, 
based on lowest order perturbation theory in the weak coupling, is independent of the PDF's 
and implies that the inelastic cross section cannot rise indefinitely. In fact, it cannot rise for 
more than about two decades in energy above the region now covered by HERA. Thus, either 
there will need to be a dramatic change in the way parton distribution functions scale not too 
far above the HERA data, or the assumption of weak perturbation theory is wrong. In either 
case the inelastic cross section will likely be very different at high Ei, than what is predicted in 
the references above. 

In Section 2 we discuss the relation between the cross section at high E^ and the low x 
behavior of the distribution functions. In Section 3 we derive our bound and Section 4 contains 
conclusions. 



II. ULTRAHIGH-ENERGY NEUTRINO-NUCLEON CROSS SECTIONS 

At ultrahigh neutrino energies the total cross section is completely dominated by its deep 
inelastic component. The deep inelastic scattering (DIS) cross section can be evaluated within 
the QCD-improved parton model, employing the universal parton distribution functions (PDF's) 
of the nucleon |T^|-p!^ which are derived from fitting photon-exchange dominated DIS measure- 
ments. Most important for high energies are the recent HERA low-x measurements [0, because 
UHE neutrino cross sections receive a dominant contribution from the ultrasmall-x region of DIS 
- even below the range covered by HERA. 

The total neutrino-nucleon cross section is given by 
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2MnE^) = dx dQ^ d^a'^^^^/dxdQ^ , (1) 
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with the standard (weak) DIS cross section of the form 



where y = Q'^/{xs). The effective quark and anti-quark distributions q{x,Q^) and q{x,Q^) 
include the appropriate (electro-) weak couplings for charged or neutral current DIS. In what 
follows we will always assume an isoscalar nucleus [A^ = {p + n)/2]. Note, that Eq. (|l|) is well 
behaved only for neutrino DIS where - contrary to the case of photon exchange - the heavy 
vector meson propagator in Eq. (H) is non-singular at = xys = 0. Accordingly, the following 
arguments are tied to the weak interaction framework, with no analogue for the more familiar 
photon exchange process. For explicit calculations we restrict ourselves to neutrino, rather than 
anti-neutrino, scattering, with the understanding that sea-quark dominance in the UHE limit 
does not discriminate between neutrinos and anti-neutrinos. 

Recent evaluations of UHE cross sections have extrapolated PDF's below x < 10~^ using 
power laws 0, radiative renormalization group behaviour or BFKL-type resummations 0. 
Within the accuracy required for astrophysical phenomenology these approaches agree and pre- 
dict a steep rise of the low-x parton density functions. Any flattening of the PDF's at low x due 
to gluon recombination is naively expected to be absent at the UHE scale Q ~ Mw due to the 
higher twist nature 0{1/Q'^) of such recombination channels [p!5| . 

To gain a feeling for the interplay between the UHE cross section and the low-x behaviour of 
the PDF's we can perform the integral under the approximation that we neglect the evolution 
of the PDF's. We findQ 



aZl is] ~ ■ / dx 
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( T— : 1 + ^(^) ( T + 1 - 2 {^-^\ ln(l + s) 



(3) 



where s = sx/M^^ and the PDF's are evaluated at some fixed scale {i.e. q{x) = q{x, M^^), 
q{x) = q{x, M^r z))- From (H) it is straightforward to show that assumed behaviour of the PDF's 



^Sub-leading QCD corrections are implicitly included in the evolution of q{x) = q{x,Q'^),q{x,Q'^); 
explicit 0{al) corrections to Eq. (||) are negligible 0(1 — 2%). 

^ i.e. imposing a "Q >few GeV" cut-off on the integral in (|I]) changes its value only marginally. 

^We have verified that the effects of evolution induce corrections of only about 20%. 
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at X — > determines the ^ oo asymptotic behaviour of the total cross section [|T6[: A 
power-type x ^ behaviour xq{x) oc 0] of the sea-quark PDF translates into a power-type 
rise with energy (y^^~*'^ oc E^. Renormalization group evolution from a fiat low-scale input 
H gives rise to a slower, but still significant, growth with energy. This is easiest to see in the 



double-leading- logarithm approximation |]T^, in which the cross section can be seen analytically 
[ITBI to rise less fast than any power but faster than a logarithm. Although each of the different 
approximations of the small-x behaviour gives different quantitative predictions of the E^ ^ oo 
asymptotics, we observe in Fig. that they agree to a good approximation up to energies as high 
as El, = 10^^ GeV. In particular, they all rise continuously with energy. 

This rise, which is a direct consequence of the x behaviour of q{x), q{x), reflects mainly 
the HERA measurements of ~ a; Z]e^(g + q){x) at small-x. The question is, therefore, up 
to what energy is the rise of crtot strictly required by these data? The ep collider HERA runs 
at a center of mass energy of ^/s ~ 300 GeV corresponding to an equivalent neutrino energy of 
^HERA ^ 5x 10^ GeV However, leaving theoretical extrapolations aside, atot for E^, > E^^^^ 
is not completely unknown because HERA probes the differential da/dxdQ'^ in Eq. (|^) down to 
much lower x (corresponding to higher E^,) than the average (log^^Qx) ~ —1.5 in the integrated 
crtot(-E'™^^), albeit at a lower scale Q < Mw,z for the 7*-exchange process. Hence, the maximal 
Ey up to which the evaluation of atot according to Eqs. (|l]) and (0) is fixed by HERA depends on 
the extent to which we trust the NLO QCD scale-evolution of DIS structure functions. Recent 



estimates |T^ of NNLO corrections to parton evolution confirm a trustworthy stability of the 
full singlet evolution at least down to x ^ 10^^ for ^ 10 GeV^. This latter kinematical 
range is well covered by HERA and the corresponding data are incorporated in the PDF sets of 
Refs. [1T2|-|14|. Assuming nothing moreQthan the universality of the PDF's and the validity of the 



RGE within 10"^ ^ x, 10 GeV^ < < M^, we can then consider (Xtot ^ if evaluated according 
to Eqs. (H) and @} - as effectively covered by HERA data up to neutrino energies where Utot 
becomes sensitive to x < 10~^. This sets inQ smoothly above Ej^ ^ 10^ GeV, which we indicate 
by the vertical line in Fig. 1. Note, this is a very conservative estimate of the impact of the 
HERA data. Were we to trust the PDF-mediated mapping of the HERA measurements onto 
the neutrino cross-section down to x ^ 3 x 10^^, and rely on the evolution from = 1.5 GeV^ 
to M^, then the most recent HERA data would cover atot at E^, ~ 10^ GeV except for 



a correction of ~ 25% from x < 3 x 10 ^. As we will deduce in the next section, this would 
suggest HERA 7-exchange data are already probing a perturbative unitarity bound in UHE 



^I.e. assumptions which are intrinsic to the perturbative evaluation in Eq. (|l|). 
^For illustration, see Fig. 3 in 0] or Fig. 3 in H. 



4 



neutrino scattering. However, for now, we prefer not to speculate about shifting the vertical 
line in Fig. 1 to the right. A detailed statistical analysis, though, might do better and provide 
stringent correlations between HERA data and Utot at even higher E^,. Note, our limitation 
X > 10~^ is also a safe condition for the absence of recombination corrections [|T5| at HERA 
scales which are, accordingly, even more negligible at the UHE scale Q = Mw,z- 



III. PERTURBATIVE UNITARITY 



Using unitarity of the S'-matrix, we can relate the total uN cross section in (|l|) to the imagi- 
nary part of the neutrino-nucleon forward scattering amplitude: 



^/A<„^-^^.) = Im 



T{uN,uN}is, t = 0) 



(4) 



where A = (s — M^)^ and s,t are the standard Mandelstam variables. The elastic amplitude T 
is related to the elastic cross section 



Jt 



1 
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T{uN,uN}{s, t) 



(5) 



Combining Eqs. (H) and (|^) gives a general limit on atot in Eq- (HD 



dt 



t=o 



16n 



0)])' 

A 



16n 



'tot 



(6) 



We note that the inequality Eq. (H), as derived, holds strictly for each spin and isospin state of 
the nucleon; however, it is straightforward to show that| 



(7) 



where the sum is over spin and isospin of the nucleon. Thus, we can just as well use the 
spin/isospin averaged cross sections in Eq. 



^The difference between tlie left and right hand sides of Eq. is the variance of elementary probability 
theory and therefore positive. 

''in fact, due to constraints on spin and isospin asymmetries in DIS scattering, one would expect the 
difference between the right-hand and left-hand sides of Eq. (|^) to be negligible in the high energy limit. 
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The inequality in Eq. (j^) is a standard statement of the optical theorem and can be found 
in many textbooks. It follows strictly from the positivity of (ReT)^ and does not rely on any 
perturbative expansion. We now consider its implications in the context of a perturbative cal- 
culation in the weak coupling. Expanding the elastic amplitude to lowest power in the weak 
coupling g and using the most general Z-nucleon coupling (see Fig. we obtain 
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4 cos 9 



w 



y M,(A;07"(l-75)w.(A;) 



t 



X Un{,p') fl{.q^) + «Cr„/3 /2( 

75 / 6-2 (g^) + qa 75 fl-sl?^ 



un{p) + 0{g') . 



(8) 



In the forward direction only the form factors /i(0) and gi{0) contribute. Using the standard 
model Z-boson current and isospin symmetry, the values of the form factors at = can 
be expressed in terms of the nucleon isospin operator T3 and charge operator Q as /i(0) = 
(2 T3 — 4 sin^ 9w Q) and gi (0) = —2 T3 g^, where gA = 1.27 is the measured value |^ of the axial 
vector coupling constant. For an isoscalar nucleus we obtain 



da. 



(2) 
cl 



dt 



t=0 



n2 



l + (l-4sin2 Owf 



+ 9I 



(9) 



Accordingly, the inequality Eq. (j^) gives 

(Ttot ^ 9.3 X 10"^^ cm^ 



(10) 



The value on the right-hand side of Eq. ([ToD is shown as a horizontal line in Fig.|l|. We observe 
that the inequality is violated for E^, ^ 2 x 10^ GeV. Note that the precise value of the bound is 
rather insensitive to the values of /i(0) and gi{0) because it depends on the square root of the 
factor in the square brackets of Eq. (^. 

HERA information on PDF's determines atot for neutrino energies as high as E,, ^ 10^ GeV. 
Perhaps large changes in the PDF's will set in within the range of x appropriate for E,y ~ 
10^ — 10^ GeV such that the bound is respected for larger E^. These corrections would presumably 
arise as a softening of the gluon and sea quark distributions at small x. 

Alternatively, one can note that the bound Eq. ( [T(]| ) is obtained by comparing an 0{g'^) 
expression for atot with an 0{g'^) expression for the forward elastic scattering amplitude. Since 
unitarity is an order-by-order statement within perturbation theory, to be strictly rigorous we 
must also include the absorptive 0{g'^) contribution to the forward scattering amplitude. In this 
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way, unitarity in the ultrahigh energy hmit can in principle be restored, but at the cost of large 
0{g^) corrections to the forward amplitude. 



For illustration, we write the elastic amplitude symbolically as 



where we have calculated the T^^^ part above and the T*^^^ term results from the exchange of a 
pair of virtual Z^^s or l^^'s in the elastic scattering, as illustrated in Fig.^. Using this notation 



da. 



ol 



dt 



oc 



t=0 



and 



^4^(2)2 ^ 2(^6^(2) Re[TW] + [(Re[T(^)])' + (lm[T(^)] 



(atot)V/ (lm[TW])' + 



(12) 



(13) 



From these equations we see that the unitarity condition Eq. (^ can be satisfied; however it 
requires that the higher order g^T^^^ term become larger than the lower order g'^T^'^^ term in 
the high energy limit. So, in this approach, calculating to lowest order in the weak coupling g 
is not a good approximation for the elastic differential cross section at high energy. But if the 
perturbation expansion is unreliable in Eq. ([T2|), then it is also reasonable to expect higher order 
corrections to ImT to be larger than the term shown, both in Eq. (p!2D and in Eq. (^3|). Since 
this is the term on which all calculations have been based it implies that these calculations may 
not be reliable at very high energies. If the perturbation expansion doesn't work in Eq. (|I2|) then 
we must question its validity in Eq. (p!3|). 



IV. CONCLUSIONS 



We compared perturbative predictions for neutrino-nucleon total cross sections at ultrahigh 
neutrino energies to a unitarity bound derived from the corresponding elastic neutrino scatter- 
ing amplitude in the forward direction. At the non-perturbative level the bound is absolutely 
stringent. A lowest order expansion in powers of the weak coupling leads to a bound which 
is saturated at an energy surprisingly close to what is effectively covered by HERA measure- 
ments. Thus, for the largest energies relevant to neutrino astronomy, sizable PDF corrections 
could reside in crtot- Alternatively, the large apparent 0{g^) total cross section should manifest 
itself through the corresponding term in the elastic amplitude as well. But if the perturbation 
expansion is not valid for the elastic amplitude how can we trust it for the inelastic cross section? 
In this case also (Xtot could be very different for large E^, than the curves shown in Fig.|l|. 
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We have ignored the possibihty that unitarity is satisfied by a new interaction - physics 
beyond the standard model p|-p!l| . But even in this alternative our basic conclusions remain the 
same: something dramatic must happen at energies close to the current energy and the existing 
calculations of atot should not be trusted at large neutrino energy. 

If the bound is correct, and the PDF's change at low x, then atot will be smaller at high 
than what is shown in Fig.[^. If perturbation theory has broken down then o"tot could be anything 
- perhaps even large enough to explain the cosmic ray data without the need for new physics. 
Or perhaps the values shown in Fig.|I| could turn out to be correct; maybe weak perturbation 
theory only breaks down for the elastic cross section. Without a better knowledge of the PDF's 
at low X, there is no way to know. 
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FIG. 1. Perturbatively evaluated QCD neutrino- nucleoli total cross sections from a power law (GQRS)[|, 
renormalization group evolution (GKR)[^, or BFKL (KMS)Q approach towards ultrasmall-x structure functions. 
The thick horizontal line is the lowest order perturbative unitarity bound in Eq. ([l0|). The vertical line is a 
conservative upper bound on the E'^-range effectively covered by HERA, i.e. the perturbative cross section for 
energies E^, < 10^ GeV receives no significant contribution from the x < 10"'* regime of DIS. 
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FIG. 2. Lowest order [0(5^)] elastic amplitude T^uN,uN} corresponding to Eq. 
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FIG. 3. Next order elastic Amplitude T^uN,uN} corresponding to T^^^ in Eq. ([iTl) . 
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